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\f*) We study the spin i hadronic state in quenched lattice QCD to search for a 

, possible S = +1 pentaquark resonance. Simulations are carried out on 8 3 X 24, 

. 10 3 X 24, 12 3 X 24 and 16 3 X 24 lattices at (3=5.7 at the quenched level with 

f*^ ' the standard plaquette gauge action and Wilson quark action. We adopt two 

I ' independent operators with 1 = and J p = | to construct a 2 X 2 correlation 

, matrix. After the diagonalization of the correlation matrix, we successfully obtain 

i the energies of the ground-state and the 1st excited-state in this channel. The 

volume dependence of the energies suggests the existence of a possible resonance 
state slightly above the NK threshold in / = and J p = ^ channel. 



1. Introduction 



After the discovery 1 of 9 + (1540), the property of the particle is one of the 
central issues in hadron physics. While it is very likely that the isospin of 
+ is zero, the spin and the parity of + and the origin of its tiny width 
of a few MeV still remain open questions. 2 In spite of many theoretical 
studies on <d + , 2 the nature of the exotic particle is still controversial. 

As for the theoretical approaches, the lattice QCD is considered as one 
of the most reliable tools for studying the properties of hadron states from 
first principle, which has been very successful in reproducing the non-exotic 
hadron mass spectra. 3 Up to now, several lattice QCD studies have been 
reported, which mainly aim to look for pentaquarks. However, the results 
are unfortunately contradictory with each other. On one hand, in refs. 4,5 , 
the authors conclude the parity of + is negative. On the other hand, in 
ref. 6 , + in the positive parity channel is reported. In refs. 7,8 , even the 
absence of 9 + is suggested. 

A difficulty in the study of + using lattice QCD comes from the con- 
taminations of the scattering states of nucleon and Kaon, since the mass 
of 9 + is above the NK threshold. In order to verify the existence of a res- 
onance state, we need to properly discriminate the first few lowest states, 
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pick up one state as a candidate, and determine whether it is a resonance 
state or not. 

In this paper, we study 1 = and J = \ channel in quenched lattice 
QCD. With the aim to extract the possible resonance state slightly above 
the NK threshold in this channel, we adopt two independent operators 
with 1 = and J = \ and diagonalize the 2x2 correlation matrix of the 
operators. After the successful separation of the states, we investigate the 
volume dependence of the energy so that we can distinguish the resonance 
state from simple scattering states. 



2. Formalism to separate the states 

In lattice QCD calculations, the energies of states can be measured using 
the operator correlations. Let us consider the correlation (O(T)O(O) 1 ') with 
0(T) an interpolating operator on the time plane t = T with a certain 
quantum number. The correlation (O(T)O(O)^) physically represents the 
situation in which the states with the same quantum number as 0(t) are 
created at t = and propagate into Euclidean time direction and finally 
are annihilated at t = T. Taking into account the fact that the created 
state vac) can be expressed by the complete set of QCD eigenstates as 
0^\vac) = cq\0) + Ci|l) + .. = $Zcj|i), we can write down the correlation 
(O(T)O(0)t) as (O(T)O(0)t) = £. £. c iCj (i\e-» T \j) = J2i \ci\ 2 e~ EiT in 
terms of the QCD transfer matrix e~ HT , with QCD Hamiltonian H and n- 
th excited-state \n) normalized as (m\n) = 5 mn . We can extract the ground- 
state energy E by taking large T so that the correlation (O(T)O(O)^) 
behaves as a single-exponential function (O(T)O(O) 1 '} ~ \c \ 2 e~ EaT . It is 
however rather difficult to obtain the excited-state energy Ei(i > 0) from 
the multi-exponential function of T, which is a problem for extracting the 
possible pentaquark signal above the threshold. To overcome this difficulty, 
we adopt the variational method using correlation matrices of independent 
operators. 9 In this method, we need a set of independent operators {O 1 }. 
We define the correlation matrix C/j as Cjj(T) = (O 1 (T)O j ^ (0)} . Then, 
from the product C _1 (T + l)C(T), we can extract the energies {Ei} as 
the logarithm of eigenvalues {e Ei } of the matrix C _1 (T + l)C(T). In the 
case when we prepare N independent operators, the correlation matrix is 
an N x N matrix and we can obtain N eigenenergies {Ei} (0 < i < N — 1). 
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3. Lattice set up 

As interpolating operators {O 1 }, we take two independent operators; 
Q 1 (x) = e abc [uJ(x)Cj 5 d b (x)}u e (x)[J^(x)j 5 d e (x)} - (u <-► d) which 
is expected to have a larger overlap with 6 + state, and 2 (:r) = 
e abc [uJ(x)Cj5db(x)]u c (x)[s^(x)j5d e (x)] — (u <-> d) which we expect to have 
larger overlaps with NK scattering states. Both of them have the quan- 
tum number of (I, J) = (0, |). Here, Dirac fields u(x), d(x) and s(x) are 
up, down and strange quark field, respectively and the Roman alphabets 
{a,b,c,e} denote color indices. In this case, the correlation matrix is a 2 x 2 
matrix. Furthermore, for a source, we construct and use "wall" operators 
®waiiW an d ©waiiW defined using spatially spread quark fields X)x9( x ) as 

©wall!*) = (v^) 5 Ex^x- 5 £ QbC ["a(^l)C754(x 2 )] Ue ( a ;3)[^(a ; 4)754(x5)] - 

(u <-> d). We fix the source operator wa ii(£) on t=6 plane to re- 
duce the effect of the Dirichlet boundary on t=0 plane, 10 and calculate 

c"(T) = E s (e / (x,T + 6)eln(6)). 

We adopt the standard plaqucttc (Wilson) gauge action at 
(3 — 5.7 and Wilson quark action with the hopping param- 
eters as (k m ,k s )=(0.1600, 0.1650), (0.1625,0.1650), (0.1650,0.1650), 
(0.1600,0.1600) and (0.1650,0.1600), which correspond to the current 
quark masses (m„ jd , m s ) ~ (240, 100), (170, 100), (100, 100), (240, 240) and 
(100, 240), respectively in the unit of MeV. We perform lattice QCD calcu- 
lations at j3 — 5.7 employing four different sizes of lattices, 8 3 x 24, 10 3 x 24, 
12 3 x 24 and 16 3 x 24 with 2900, 2900, 1950 and 950 gauge configurations. 
At (3 = 5.7, the lattice spacing a is set to be 0.17 fm so that the p meson 
mass is properly reproduced. Then, in the physical unit, the lattice sizes 
are 1.4 3 x 4.0 fm 4 , 1.7 3 x 4.0 fm 4 , 2.0 3 x 4.0 fm 4 and 2.7 3 x 4.0 fm 4 . 

We take periodic boundary conditions in the directions for the gauge 
field, whereas we impose periodic boundary conditions on the spatial di- 
rections and the Dirichlet boundary condition on the temporal direction 
for the quark field. Many lattice studies about + adopted a periodic 
or antipcriodic boundary condition on the temporal direction for quarks. 
In the case when the temporal length N t is not long enough, we have 
to be careful about the contaminations by particles which go beyond the 
temporal boundary. For example, with the periodic boundary condition, 
the correlation (0(T)0(O)) contains the unwanted contributions such as 
(K\Q(T)\N)(N\Q(0)\K) ~ e -E N T+E K (T-N t ) _ whilc n0 quark can g0 over 

the boundary on t=0 in the temporal direction, the boundary is transpar- 
ent for the particles composed only by gluons; i.e. glueballs, due to the 
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periodicity of the gauge action. Taking into account that these particles 
are rather heavy, it would be however safe to neglect these gluonic parti- 
cles going beyond the boundary. Then, the correlation (0(T)0(O)) mainly 
contains such terms as (vac\Q(T)\N K) (N K\Q(0)\vac) and we can use the 
exactly same prescription mentioned in the previous section. 



4. Ground-state and 1st Excited-state in (J, J p ) — (0, | ) 
channel 



In this section, we investigate the ground-state and the 1st excited-state 
in (7, J p ) = (0, i ) channel. We focus on the volume dependence of the 
energy of each state, to distinguish a possible resonance state from NK 
scattering states. We expect a resonance state to have almost no volume 
dependence, while the energies of NK scattering states are expected to scale 
as yj Mjj + |^n| 2 + yj M|- + |^n| 2 according to the relative momentum 
^n between N and K on finite periodic lattices on the assumption that the 
NK interaction is negligible. Note here that there are other possible finite 
volume effects due to the finite size of N and K or possible pentaquark 
states, than the volume dependence arising from ri. We then need to 
take account of this fact in the following analysis. 

Now we show the lattice QCD results of the ground-state and the 1st 
excited-state in I = and J p = | channel. In Fig. 1, we plot the ground- 
( K u,d. K s)=(°- 1 650,0.1 650) (K ud ,K s )=(0.1 600,0.1 600) 
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Figure 1. The lattice QCD data of the ground-state in (I,J P ) = (0, ^ ) channel arc 
plotted against the lattice extent L. The left figure is for (k u d> K s) = (0.1650,0.1650) 
and the right figure is for (k Uj£ j, k s ) = (0.1600, 0.1600). The solid line denotes the simple 
sum Mjv + Mx of the masses of nucleon Mpf and Kaon Mji . We use the central values 
of Mjv and Mj[ obtained on the lattice with L = 16 
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state energies in I = and J F = \ channel on four different volumes. 
Here the horizontal axis denotes the lattice extent L in the physical unit 
and the vertical axis is the energy of this state. The symbols with error-bars 
are lattice data and the solid line denotes the sum Mjv + Mk of the nucleon 
mass Mm and Kaon mass Mk- At a glance, we find that the energy of this 
state takes almost constant value against the volume variation and coincides 
with the simple sum Mjy + Mk- We can then conclude the ground-state 
in I = and J p = | channel is the NK scattering state with the relative 
momentum |p| = 0. 

(K ud ,K s )=(0.1 650,0.1 650) (K ud ,K s )=(0.1 600,0.1 600) 




Extent L [fm] Extent L [fm] 

Figure 2. The lattice QCD data of the 1st excited-state in (I, J p ) = (0, ^ ) channel are 
plotted against the lattice extent L. The solid line represents y + pp+y'' M^- + |p| 2 
with |p| = 2tt/L the smallest relative momentum on the lattice. The dashed line denotes 
M N + M K . 

Fig. 2 shows the 1st excited-state energies in this channel again in terms 
of L. The solid line shows y/ M N + {2n/L) 2 + \J M K + (2tt/L) 2 as the 
function of L. In this case, the lattice data exhibit clearly some volume 
dependence. One would naively consider this dependence to be that of 
the 2nd lowest NK scattering state, which scales as \JM N + (2ir/L) 2 + 
\J M K + (2tt/L) 2 , and would take the deviations as from other possible 
finite size effects. However, if it is the case, the significant deviations in 
1.5 < L < 3 fm in the right figure may not be justified. In the range of 
L where the lattice data in the left figure scale just as the solid line, wc 
should observe the same behavior in the right figure. Because, in the case 
when quarks are heavy, composite particles will be rather compact and we 
expect smaller finite size effects. 

We can understand this behavior on the assumption that this state is 
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a resonance state rather than a scattering state. In fact, while the data in 
the left figure rapidly go above the solid line as L is decreased, which can 
be considered to arise due to the finite size of a resonance state, the lattice 
data exhibit almost no volume dependence in the right figure especially in 
1-5 ^$ L < 3 fin, which can be regarded as the characteristic of resonance 
states. 



5. (I,J P ) = (0, i + ) channel 
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Figure 3. The lattice QCD data in the (I,J P ) = (0, \ + ) channel are plotted against 
the lattice extent L. The solid line denotes the simple sum Mjv* + of the masses of 
the ground-state negative-parity nucleon Mjv* and Kaon Mk- 



In the same way as (/, J p ) = (0, \ ) channel, we have attempted to 
diagonalize the correlation matrix on (7, J p ) = (0, ^ + ) channel using the 
wall-sources Owaii (i) • In this channel, the diagonalization is rather unstable 
and we find only one state except for tiny contributions of possible other 
states. We plot the lattice data in Fig. 3. One finds that they have almost 
no volume dependence and that they coincide with the solid line which 
represents the simple sum Mm* + Mr of M^* and M^, with M/v* the 
mass of the ground-state of the negative-parity nucleon. From this fact, 
the state we observe is concluded to be the N*K scattering state with the 
relative momentum |p| = 0. It is rather strange because the p-wave state 
of N and K with the relative momentum |p| = 2n/L will be lighter than 
the N*K scattering state with the relative momentum |p| = 0. We miss 
this lighter state in our analysis. This failure would be due to the wall- 
like operator 6 W aii(£) : The operator 6 W aii(£) is constructed by the spatially 
spread quark fields X)x9( x ) with zero momentum. This may lead to the 
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large overlaps with the scattering state with zero relative momentum. It is 
also desired to try another operator which couples to p-wave NK scattering 
state. 

6. Summary 

We have performed the lattice QCD study of the S=+l pcntaquark baryons 
on 8 3 x 24, 10 3 x 24, 12 3 x 24 and 16 3 x 24 lattices at 0=5.7 at the quenched 
level with the standard plaquette gauge action and Wilson quark action. 
With the aim to separate states, we have adopted two independent opera- 
tors with I — and J p = | so that we can construct a 2 x 2 correlation 
matrix. From the correlation matrix of the operators, we have success- 
fully obtained the energies of the ground-state and the 1st excited-state in 
(I, J p ) = (0, \ ) channel. The volume dependence of the energies shows 
that the 1st excited-state in this channel is rather different from a sim- 
ple NK scattering state and is likely to be the resonance state located 
slightly above the NK threshold, and also indicates that the ground-state 
is the NK scattering state with the relative momentum |p| = 0. As for the 
(I, J p ) = (0, \ ) channel, we have observed only one state in the present 
analysis, which is likely to be a N*K scattering state of the negative-parity 
nucleon N* and Kaon with the relative momentum |p| = 0. 
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